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Abstract. A unique form of turbulent-transport equations is derived based on first 
principles. The role of noncquilibrium statistical mechanics employed to describe the 
phenomenology is that it enables to single out the unique form consistent with master 
equation of Liouville, a prerequisite not met with existing equations for turbulence 
modeling. The equation is variable-separated to yield a Navier-Stokes equation 
in 6D (physical-plus-eddy) space with homogeneous boundary conditions. Turbulent 
transports such as Reynolds' stress are calculated using a solution of this equation; 
a solitary-wave function. Satisfactory agreement is observed with existing experiment 
for mixing shear layer of incompressible flows although no empirical constants to fit 
with data are involved. 



1. Introduction 

First principle bases of phenomenologies of fluid dynamics and thermodynamics are 
due, respectively, to Chapman- Enskog|]lJ|,pl and Prigogine|§, who showed the validity of 
Navier-Stokes equation and the equality expression of the second law of thermodynamics 
on common basis of the Boltzmann equation expanded to the first order deviation from 
equilibrium. 

It was, however, an open question whether these equations still hold for 
turbulent flows where the dependent variables are stochastic and fractal, therefore not 
differentiable||]. This problem was solved[|] using the microscopic density |J, namely, 
unaveraged Boltzmann function, leading to rederivation of the Navier-Stokes equation 
written in terms of instantaneous quantities without assuming any statistical concepts 
like local equilibrium, or first order deviation form it. Thus the equations currently 
employed for the direct numerical simulation(DNS) has acquired the first principle basis. 
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Then a question will arise as to whether the alternative methodology of 
computational fluid dynamics, namely, Reynolds- average formalism can be founded on 
first principles as well. 

In contrast with DNS founded on the microscopic density(the Klimontovich 
formalism) its Reynolds-average counterpart should be based on its averaged version, 
namely, the Boltzmann function. Structure of this Boltzmann formalism, so to say, is 
shown to be identical with what is called the BBGKY hierarchy theory of nonequilibrium 
statistical mechanics 

Each of the two formalisms has advantages as well as disadvantages from 
computational viewpoints. The most serious disadvantage of the K-formalism lies in 
the fact that it deals with fractal quantities which are selfsimilarly rugged to a length 
scale as small as the Kolmogorov scale ||. It means that if the 3-D Navier-Stokes 
equation is to be solved by a finite difference method, the grid size be smaller than this 
length, namely of 0(R~ 3 ^) (R ; Reynolds number). It requires the computer memory 
size growing with i? 9 / 4 (the small eddy difficulty). For this reason, the consensus upper 
limit of applicability of the current DNS falls short of R ~ 10 4 . To be able to apply 
the DNS for practical design of transport vehicles(i? ~ 10 7 ) we would have to have a 
computer with memory size greater by the factor of (10 7 /10 4 ) 9 / 4 ~ 10 7 compared to the 
currently available ones. This situation is not changed even when one employs spectral 
methods to avoid the conceptual difficulty of having to work with differential equations. 
For, then, the number of the Fourier modes to be taken into account increases at the 
same rate. 

The B-formalism, in contrast, is free from such difficulty owing to the statistical 
average taken in the process of generating distribution functions at the expense of 
dealing with multitude of such functions, infinite in number. The statistical average 
as meant here is either an ensemble average over repeated experiments, or an average 
over time that is long enough for the fractal ruggedness to be smoothed out, yet is short 
enough for fluid-dynamic unsteadiness, such as shedding period of Karman vortices or 
aerodynamic flutter to be discernible. The ergodic theorem warrants their identity. 

The objective of this paper is to demonstrate that the Reynolds average regime 
of the computational fluid dynamics can be founded on a firm first principle ground 
through the following processes : 

i) Demonstration of how a closed set of kinetic equations is obtained out of an infinite 
chain of them of the B-formalism. 

ii) Taking fluid moments of the kinetic equation for two-point fluctuation-correlations 
that is unique from the viewpoint of first principles. 

iii) Rederivation of the whole set of equations thus obtained from phenomenological 
equations being used for DNS, namely, eliminating the process via statistical 
mechanics. 



2. The Boltzmann and Klimontovich formalisms : A review 



3 



Any statistical theory rests on the axiom that a field quantity f(z) describing stochastic 
and possibly fractal physical phenomenon is equivalent to a set of quantities that are 
smooth, deterministic, and infinite in number, 



/(=/) 



/'/' 



f'f'f 



(/' = /- /) 



(1) 



where overbar stands for the statistical average as defined in the preceding section, and 
/ = f(z) is the same quantity at a different point z. 

The Klimontovich formalism as defined here is a formalism where / is identified 
with the microscopic density || 

l(z)= £ 8{z-zW(t)) (2) 

l<n<AT 

In this expression z = (a;, v) is a point in the phase(/i-) space, z^ n \t) is the locus of n-th 
molecule in this space, N is the total number of molecules under consideration, and 8 
denotes the six-dimensional delta function. 

It has been shown || that the equation governing / is the unaveraged Boltzmann 
equation; 

d 8 



(3) 



with J denoting the classical collision integral acting on molecule z . The key issue here 
is that Eq.(|3]) is a deterministic equation of continuity in the /z-space, free from any 
statistical concept s|5]]. 

The Boltzmann formalism, on the other hand, deals with quantities to have appeared 
on the r.h.s. of (1), subject to statistical average, namely, the Boltzmann function 



and the correlation functions of consecutive hierarchies 

ip(z, z) = 



(4) 



f'f 



ip{z,z,z) = f'f'f 



f f-ff 



f f f-Mz,z)-Mz,z)-Mz,z) 



(5) 
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Klimontovich variables || such as instantaneous gas density p and center-of-mass 
velocity u of the iV-particle system are given by 



{■OO 

p = m 



/oo 
/ dv = m ^ 6(x — x^ n '(t)) (m : mass of a molecule) (6) 
-OO , x „ x >r 



Kn<N 



/oo 
vf_dv = m v (n \t) 6{x-x (n) {t)) (7) 

-OO ix_^Kr 



Kn<N 



The Reynolds-averaged version of those fluid variables are generated by taking average 
of expressions (f|) and (0) 

/oo 
/ dv (8) 
-oo 

/oo 
vf dv (9) 
-oo 

Two-point correlation i/j(z, z) consists of two parts; short-range part due to direct 
molecular collisions that is irrelevant to turbulence and long-range part attributable to 
turbulence correlations. The latter part is expanded in a double series of Hermite 
polynomials ||10|| where all the turbulent transport terms including Reynolds' stress 
appear as expansion coefficients. 

A set of equations governing those quantities standing on the l.h.s. of (|5]) is 
generated from Eq.(|3|) by taking moments and averaging. This procedure is not 
unique; in general there are infinite ways of constructing such moment equations. This 
arbitrariness is eliminated by invoking a postulate that the whole system be consistent 
with Liouville 's equation(the equation of continuity of iV-particle probability density in 
6iV-dimensional space; the master equation that is universally valid), or its corollary that 
the whole set of equations be identical with those of the BBGKY theory at each level of 
hierarchies. The only difference is that the BBGKY generates destribution functions in 
the direction of descending number of molecules through a series of integrations starting 
from iV ~ O(10 20 ), whereas here are defined the same functions in ascending number 
of molecules. The kinematical information missing in the latter approach is identical 
in analogy with the fact that we cannot predict the functional form of g{x,y) out of 

/oo 
g(x,y)dy . This is why we need the postulate, 
-oo 

The averaging process consistent with the postulate has led to following set of equa- 
tions H : 



1 - particle level 



d d 



fi [/] = 0: (^ + «-^)/-W)[// + #*,2)]=0 (10) 
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2 - particle level 



f'B[f) + f'B[l)=0 , (B[ ] = {B[ 

^(2, z, z) ] - J (2(2) [ /V>(*, 2) + /y>(*, z) + 2, z) } = 
3 - particle level 



(11) 



/7'^[/] + ff'B[f_\ + /'/'£[/] = (12) 

This system constitutes a chain of equations for the set of infinite number of variables 
[f(z),ip(z,z),i/;(z,z, z), ■ ■ •]. 

The issue that is most crucial to the quality of the proposed approach is the closure 
condition. It is how to truncate the infinite chain of equations to make the system 
tractable without violating physical soundness. 

Early stage of development along this line has employed the following condition of 
tertiary chaos [|T(J, fl~T| 



ip(z,z,z) = (13) 

a condition next to the simplest one known as Boltzmann's (binary) molecular chaos 
hypothesis; i/;(z,z) = 0. 

Eq.flTTD under this closure condition is investigeted in some depth: It is shown 
that assumption (|IJJ) allows Eq.flTTD f° r separating variables into those for respective 
particles, thereby its fluid moment equation leads to linearized Navier-Stokes equation 
(the Orr-Sommerfeld equation, in particular.). It is also seen that this closure gives 
satisfactory description only for weak turbulence where the nonlinearity in turbulent 
intensity does not play major roles. Later, it has been superseded by alternative one 
that has wider range of applicability, yet preserving the variable-separability beyond 
linear regime ||12j| : Namely, put 



30 



tp(z,z)—K.P.r (p(z,(v)(j)*(z,uj) duj 
J — 00 

/oo 
(j)(z, lu)4>(z, Q)5(lu + 0) dujduj (14) 
-00 



ip(z, z, z) = R.P. t 2 / <f>z{z, u)4>s(z, Q)(f) 3 (z, uj)5(lu + uj + uj) dujduduj (15) 
J — 00 

where uj is the variable-separation parameter having the dimension of the frequency, r 
is a characteristic time, symbols(*) and R.P. denote the complex conjugate and the real 
part, respectively. Separated variable is complex, subject to 

<f>*(z,u) = <f>{z,-u) (16) 
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as will be justified a posteriori. (See Eg. (|T8|) below.) The closure condition is introduced 
in the following form 

03 = (17) 

It makes Eq. (|TTD separated into two equations each for respective points z and z, in the 
form of complex conjugate to each other provided that condition (|T6|) is met : 

(-iu + — + v • — ) - J(z|s) [0/ + f4> + r / 0(w - uj)4>(uj)duj] = 0(18) 



<9t 9a;' 

If the nonlinear term of convolutional integral is deleted, the equation degenerates to 
the previous case with tertiary chaos closure. 

Fluid moments of (ff^) together with (14) provide equations for turbulent transports 
such as Reynolds' stress and turbulent heat flux density to obey. For actual derivation 



of these equations see [13 



3. Fluctuation equations in physical-plus-eddy space 



The approach described in the previous section has shed some lights in turbulence 
research. In fact, for the cases tested agreement with experiments is satisfactory 
although the theory is free from any adjustable parameters M, |15l as contrast with 



existing models such as eddy-viscosity model. The success, however, has been limited 
to cases where the flow geometry is governed by single variable. (Note that the velocity 
fluctuations are multi-dimensional.) 

We will show that a small renovation of the theory sketched in the preceeding section 
can make turbulence with general three-dimensional geometry tractable. It is to replace 
the frequency u, a scaler quantity having appeared in Eqs.(|i~8|) and flOf), by the wave 
number k that is a vector connected to the frequency through phase velocity c by the 
dispersion formula 



UJ 



c ■ k 



Then, new separation rule to replace (14) is 



ip{z,z) 



R.P. /' 



(z, k) 4>*(z, k) dk 



(19) 



(20) 



where I is the characteristic length of the macroscopic phenomenon under consideration. 
Accordingly the governing equation for <p to replace (|T8D is written as 



iuxj) = Q((f>) 
with 



+ 6f + R.P. V 



[k-k)<p(k)dk 



(21) 
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For practical purposes it is more convenient to separate out the spacially periodic 
factor from by putting 

(j)(z,k) = e lkx $(z,k) (22) 

and work with its amplitude <3>. Owing to the fortuitous situation that the only nonlinear 
term in Eq.(^) has the form of convolutional integral, the factor e tk ' x is seen to drop off 
upon substituting expression (p2[) into Eq. (|21|) . Since (d/dxj)4> = e lkx (d/dxj + ikj)$ 
the equation for $ should preserve the form of Eq.(^Tj) with the only substitution 

<} <} ikj (23) 



dxj dxj 

Another favorable property of the nonlinear integral is that it reduces to a simple 
product through Fourier transform 

$(z, k) = — ^ /_~ e-* k - s F(z, s) ds (24) 

such that 

POCi 

$(fe - k) $(fe) dk — > FF (25) 
Then Eq. ([H]), with (|22| ) through (|25| ) taken into account, reads 

{lft ~~ ° 3 lk- + Vj9j ) F ~ J{Z ^ + * F + FF } =0 ^ 
where dj is six- dimensional operator defined by 



The Fourier variable s introduced through the transform ( p24|) has the dimension 
of length characterizing the size of eddies. Therefore s might well be called the eddy 
variable. 

Substituting ( p2|) and (p4|) into (|2~0D we have a remarkably simple formula for ip 
when written in terms of F, 



1 Z" 00 

Vv, z) = j2rfj3 J F(x,s,v)F(x,s + x-x,v)ds 



(28) 



Note that here all the quantities of the integrand are real in contrast with those of (|20|). 

Derivation of fluid equations from Eq.(p6[) is almost parallel to its low dimensional 
predecessor |T3[, so will not be repeated here. Minimum necessary steps to reach the final 
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form of fluid equations will be given in order : First we employ the same fluid-moment 
expansion for F(x, s, v, t) in terms of Hermite's polynomial 7i in the -u-space after Grad 

0; 



-e/2 



(27ra 2 ) 3 /2 m 



„ j_ ft -U I -L/ I Ijkk n. 



■jrr 



with 
H 

t 

u) 
rrij 



w/a = (v — u^)/a 
rrij/p 



puj + p'u'j 



(29) 



a 2 = RmT (Rm = Avogadro no. /molecular no.) 

Note that this expansion differs from the classical 13-moment expansion in that 
qj 7^ 0,qjj 7^ 0. Second, we replace Eq.(^) by its moment equivalents 



/oo 
M a [Eq.(^)]dv = 
-oo 



(30) 

where F is substituted by fl29|) and M a stands for one of the 13 moment functions; 
l,T-Cj,Hjj,TCjk,'Hjkk- These choices of moments correspond to equations of fluctuations 
of continuity, momentum, energy, stress tensor qj k and heat flux density vector qj kk , 
respectively. The actual form of the equations are : 

Dq + d r q r = 



Dqj + d r q jr + djq 40 + 



du\ 
dx 



1 dp 
Qr tj qo 



p dx 



3 5 dv) 

-Dq i0 + d r q rjj + -d r (a 2 q r ) + T^f<?40 + 



= 

du\ 

dx, 



-Qjr 



1 dp 

p dx r 



q r = o 



<?40 
qj r 



q 4 + a z q 



-p 



qk 



p 



1 dp 



q r jj 



pR M dT 
A 



' V P / 

du\ dvl 



dx T 



2 d4 

- °jr 



R 



M 



.P , 



dxj 3 dx k 
1 dX dT 



P I 

94 



1 

+ - 

P 



T^jrqk 



pRhi dT dx, 



5q r q 4 
2 p 



where we have defined the following symbols : 

Dq ^Yt' Cr d7 r +dMq) 

In the above, p, p, T, p and A denote density, pressure, temperature, viscosity and 
thermal conductivity coefficients, respectively. 



(31) 
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It is readily seen that Eqs.flaTD are exactly the same as those derived in 
ref. [|T3| [Eqs. (|21|') through (plf)] if the following replacements are effected, 



™ — <*A (32) 



d d d , . 

+ or 33 



dxj dxj dsj 

The set of Eqs.(|3"ID describes evolution of five quantities [qo,Qj,Qi) subject to 
homogeneous boundary conditions that all g's vanish with |s| — > oo, on the solid surface, 
and wherever turbulent intensity is zero in the physical space. Therefore, the expected 
solution for those quantities must have the form of a solitary wave (not to be confused 
with a soliton). 

Eqs.(|3"T|) are equations governing compressible turbulence. For incompressible flows 
(<Zo = 0) the energy fluctuation equation is decoupled, and a closed set of equations that 
results is 

d jqj = (34) 

011 ' Q 

(D - ud 2 r ) qj + d jq40 + -r-^qr + -d rgj = (35) 

ox r p 

where v is the kinematic viscosity, dj and D have been defined by (p7] ) and (|31|), 
respectively. It is readily seen that these equations represent the equation of continuity 
and the Navier-Stokes equation generalized to 6D space (x, s). In fact, if we suppress 
the eddy variables(<9/<9sj = 0) this set degenerates to the usual Navier-Stokes equation 
for velocity Uj + p~ x q^ and pressure p + q i0 . If further, nonlinear terms are neglected in 



Eq.(p5|) and parallel flow(uj = 8jiu(x2)) is assumed, they reduce to the Orr-Sommerfeld 



equation to govern p 1 q 2 as it should. 

4. Correspondence rule; relationships between solitary-wave functions and 
observables 

It should be remarked that variables q's do not correspond to any turbulent fluctuations 
that are tangible to macroscopic sensors. They are shown to be related to fluctuation 
correlations of turbulent quantities through the following deduction: From (§), (@) and 
(1) we have expression for instantaneous density fluctuation that is an observable; 

roo , 

p = ml f dv (36) 
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Similarly, by subtracting from (0), we have 



pu'j + p'uj + (p'u'j — p'u'j) =m Vjf'dv 

J—oo 

This expression, upon substituting (KB) for p', reduces to 



/oo 
Wj f'dv + 0(f 3 ) (37) 
-oo 



puj = m 
with 

0(/' 3 ) = p -VaJ-(pX-pS) 

from which we have, utilizing definition @ for ip and neglecting terms of 0(f' 4 ) 



ppu'j u[ = m r I Wj wi ip(z, z) dvdv (38) 



Furthermore, by substituting (p5) for ip and then fl29|) for F, and carrying out integration 
over(i>, v) with orthonormal property of the Hermite polynomials incorporated, the 
following relationship results; 

1 fo° 

ppu'j u\ = j^fjs J q ^ Xj s ) qi(x,s + x- x) ds (39) 
Turbulent intensity or Reynolds' stress is then given by putting x = x\ 



1 f°° 

u ' j u ' 1 = P 2( 27r /)3 y_ 00 qj ( x ' ^ ^ ds ^ 

a relationship expressing the observable turbulence intensities by an integral operation 
of the wave function. It is through this relationship that the Reynolds averaged Navier- 
Stokes equation is coupled with Eqs.(|3"ID that govern q's standing on the r.h.s. of (BUI). 



In a similar fashion the wave function representing temperature fluctuation T can be 
derived as follows: Since the ideal gas law p = RmpT_ holds for instantaneous variables 
||, we have 

R M pT' = p' - R M Tp' 

m f°° -if 1 a 2 r f >, 

= — / Wjjdv — am j dv 

3 J — OO J—oo 



ma 



— oo 
2 i-oo 



r°° , 
/ Hjjf dv, 

J — oo 



Accordingly, from fl28|) 
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where we have employed Hermite expansion Thus we see that g 4 = (1/3)^ is the 
wave function to be responsible for the temperature fluctuation R M pT' . 

Summarizing, the following correspondence rule holds between untangible wave 
function q a and corresponding observable fluctuation A' a ; 



( 



\ 



3 9jj 



V 



<?40 



5 A a 



} 



pu> 
pR M T' 



(41) 



V p' J 



where the fourth quantity q^o in the column of q a is linearly dependent on q and g4.(See 
Eq.(pT|).) They are related to each other through the following fluctuation correlation 
formula ; 

1 



(2irlf 



q a (x, s) qp(x, s + x — x)ds 



(42) 



In particular, for example, turbulent heat flux density is given as follows, using (|41|) and 
(|42T) with x = x\ 



2 rn' I 

c p p T u'j 



5 /■ v 
2 



g4 qjds 



(43) 



where c p is the specific heat under constant pressure. 



5. Reformulation within phenomenologies 



Once the correspondence rule ([IT]) has been established we are able to reconstruct 
Eq.(PTD using phenomenologies alone. This is what is anticipated because the present 
(Boltzmann) formalism is the averaged version of the Klimontovich formalism describing 
A' a directly, where the identity of phenomenologies and first-principle approach is 
warranted p. 

To effect this we shall base on the assertion that turbulent compressible flow of inert 
gas is governed by the following set of equations: 



dm, r 



dp. 
A = -= l 

dt ' dx r 







(44) 



A; 



dm a 



+ 



d 



dt dx 



ULjEkr 



+ 



E 5 jr + (E jr ) 



NS 



(45) 



d 



A =IL(e + ^ 
4 dt \~ 2p J ' dx r 



_p\-2p) p \^r) 



NS 



Q , 

.—r/ F 



0(46) 
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where 



— 3 = P—j 



(47) 



E_ = pe = p/ (7 — 1) (7; specific heats ratio) 
# = ^ = p[ 7 /( T -l)] 



(48) 



■jrj NS 



P 



d I mj 



d ( m r 



dx r \ p J dxj \ p 
A d ( p\ 



2 d (m k 

3 jr 9o; fc 1 p 



•r/ J7 



(49) 



denote, respectively, the mass flux density, the internal energy and the enthalpy per unit 
of volume (the lower case letter refers to specific quantities), the Navier-Stokes law for 
stress tensor and the Fourier's law for the heat flux density. 

It has been shown that plain (Reynolds) average of this set of equations 



A. 



(« = 0,j,4) 



(50) 



is in exact coincidence with the first principle deduction for monatomic gases (7 — |jf!7 
Written explicitly, they are 



A ^ + ^ = 

Ot OXr 



(51) 



A, 

Pj r 



drrij 



+ 



d ( mjm r 



dt dx 



P 



+ pSj r + p jr = 



(Pjr) NS + PUjU'r 



?-3 r J NS 



(52) 



_ d / mf\ d (m r TT 1 9 m,- ^\ 



Qr = (Qr) F + PCp T'u' r 



•(53) 



(Qr) F = [Q rJl 

A few remarks are in order : 

i) These equations are exact to 0(A' 2 ). 
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ii) They are written in terms of quantities that are proportional to the density, for 
example, mean mass-flux density 



rrin = m 



TUj = pUj + p'Uj 



to replace mean fluid velocity Uj, also mean internal energy E per unit of volume 

E = pe = —— (p + R M JT) 
- 7-1 \ / 



7 

to replace the specific internal energy e. 

iii) Item ii) is the key that enables to express Reynolds- averaged equations fl52|)and (j5"3|) 
in compressibility invariant forms, in other words, single term turbulence correction 
to each of stress pj r and heat flux density Q r suffices even under presence of density 
fluctuation such as turbulent combustion . 

iv) Otherwise, lengthy additional terms for turbulence correction would appear, or we 
would need the so-called mass averaging (say Uj, for instance) that suffers from a 
conceptual difficulty (u'j ^ 0^ in processing experimental data. 

Next step, the main issue of this section, is to show the identity of Eqs. ([U]) with 
the following phenomenological equations 



A' a A f3 + A' fB A a = 0, (a,/3;l,j,4) 



(54) 



where X a is defined in (|4~l"l) . These equations consist of terms of double (0(A' 2 )) and 
triple (0(A' 3 )) correlations for which we employ the separation rule exactly parallel to 
those of the previous section [Eq . (p0[) through Eq.( 
Put 



A' a A' = R.P. Z 3 / g a (x, k) g p (x, k)* dk 



oo 



A' a A'A' = R.P. I 6 / g a (x, k) g p (x, k') g 7 (x, k - k')*dkdk' 



oo 



/oo 
g a (x, k') gp{x, k)* (? 7 (5;, k — k') dkdk! 
-oo 

with g a (x, -k) = g a (x, k)* 

and substitute into Eq.(|^), then we are led to the equation that allows for the separation 
of variables : 



(55) 



dk g a (g p )* 



At 



+ 



At 



_J_ 

9(3 



(56) 



%0J — 100 
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In the above is the fluctuating part of A a in which A' a is replaced with g a , also A a A 

/oo 
g a (k') g 7 (fc — k')dk\ and d/dx r with d/dx r + ik r . The 
-oo 

separated equation thus obtained, namely, 



iujg c 



+ Aj, = (57) 



is then rewritten in Fourier-analyzed form after 

j roc » 

9{X ' ^ = j-oc e q ( X > S) dS (58) 

to lead to the same equations as Eqs. (|3TD except for the energy fluctuation equation 
which takes the form, 

1 7 dv) dxi}- 1 Op 

-Dq m + d r q rjj H -d r (a 2 q r ) + ^-^o + ^T'Qjr tt~ Qr = 0(59) 

7 — 1 7 — 1 ox r ox r p ox r 

This equation reduces in the case of monatomic gases (7 = 5/3) to the third of Eqs.(|3T|) 
as it should. 



6. A solitary-wave solution for mixing layer turbulence 

A preliminary computation checking whether the present approach is physically sound 
has been carried out for turbulent mixing shear layer of an incompressible flow [fTsTl . 
Eqs.(|34]) and (|35|) are employed assuming the average flow profile [u(r]), v (77), 0] with 
7] = X2/0JX1 as prescribed. The flow is self-similar in this sense as confirmed by 
experiment |19| , which is an indicative of molecular viscosity playing no roles in the 



equation. (See ref. [0.) Then we have the following set of equations : 

(60) 



diqi + d 2 q 2 + dq 3 /ds 3 = 
NL q 1 + d\q A Q - a-qu'qi + u'q 2 = 
NL_q 2 + d 2 q m - arjv'qx + v'q 2 = 
NLq 3 + dq m /ds 3 = 



where u' = du/drj, and 

d\ = (1 — asi)d/dsi — a(r]d/drj + s 2 d/ds 2 + s 3 d/ds 3 ) 
d 2 = d/d7] + d/ds 2 } (61) 

NL = d/dt-cd/ds 1 + (u + qi)d 1 + (v + q 2 )d 2 + q 3 d/ds 3 

with s redefined using the shear mixing layer thickness I = axi, namely, s/l — > s. 

The set of equations has five independent variables (s,7],t), therefore no existing 
tools are immediately available. At this preliminary stage of checking physical soundness 
of the proposed approach it is advisable to suppress variable s 2 by assuming d/dr] » 
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(62) 



d/ds2- The set of equations is solved for arbitrary chosen initial values for q a and 
boundary conditions 

qj — > , (J = 1, 2, 3) as |s|, \t)\ — > oo 
<?40 — > as |s| — > oo , dq^o/drj — > as 77 — ► 00 

The form of solution to be expected from these homogeneous boundary conditions 
must be a solitary wave generated by the shearing motion and kept sustained by 
nonlinearity. 

In Figs.l are shown such standing waves that build up and reach steady state with 
elapse of time for different choices of the wave speeds ; (a) experimentally observed one 
c = u = [u(oo) + u(— oo)]/2, and (b) Taylor's hypothesis c = u(r]). Reynolds' stress 
[Eq.(f40Dl is calculated using this solution and compared with existing experiment flHJ 



in Fig. 2. Agreement is more than reasonable considering that the theory involves no 
empirical constants to fit with experiments. 

7. Comparison with classical statistical theory and current turbulence 
models 

Consistency with classical statistical theory 



Karman and Howarth |p0| , the founders of classical statistical theory of turbulence, 
have derived Eq. (|54l) for incompressible flows (a, (3; 0, 1, 2, 3) correctly on intuitive basis, 
with no reference to first principles. Obviously these 6-D equations are not tractable in 
this form, so homogeneity/isotropy assumptions have been introduced. Here we have 
employed the method of separation of variables, thereby the classical limitation on flow 
geometry is eliminated. It is effected, however, at the expense of introducing additional 
independent variable s in the Navier-Stokes equation as we have seen through Eqs.(|34]) 
and (H) . 

The classical assumptions of homogeneity and isotropy are often referred to as 
an oversimplification of reality. In fact, also here, this model is shown to lead to an 
unphysical solution : For the equation of continuity 

£-° < 63 > 

which is the homogeneous version (d/dxj = 0) of Eq.(^), coupled with the isotropy 
requirement (Robertson's theorem pTf ) 

Qj = s J ( l( s ) > ( s = \ s \) ( 64 ) 



gives 



3g + s ^ = (65) 
as 
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Solution of this equation, namely, q ~ s -3 causes the integral ( JHf ) for turbulence 
intensity to diverge. To be noted is the fact that this is a direct consequence 
of the equation of continuity, a kinematical relationship universally valid, therefore 
independent of any closure condition employed. 

Inconsistency with k — e models 

We have seen that formulation in 6-D space (x, x) [Eq.(|5lD1 is the only one that 
is consistent with Liouville's equation, namely, with first principles of nonequilibrium 
statistical mechanics. Turbulent transports such as fliOD and ( pE5| ) that are quantities in 
3-D space are obtained from the solution in the 6-D space by putting x = x after the 
equations have been solved. Suppressing variables in the equations as is often employed 
is a pathological process. The following simple example would help extract what is 
meant by this trivial-looking warning : Let a steady-state temperature distribution of 
a three-dimensional body, say, a column with rectangular cross-section in the x-y plane 
be asked, but the information only on the diagonal plane x = y is required. Needless to 
say that proper process is to work with 3-D Laplace equation for the solution T(x, y, z), 
and put x — y to have T(x, x, z). The improper process mentioned here corresponds to 
solving pathological equation {2d 2 /dx 2 + d 2 /dz 2 )T = 0. 

Majority of turbulence models currently prevailing do not concur, in the simplest 
case of the homogeneous and isotropic turbulence, with the Karman-Howarth theory as 
a consequence of the hasty reduction in independent variables as sketched here. 

Summarizing, the proposed approach shares the common basis with k — e models 
only at the lowest level of description, namely, the Reynolds-averaged Navier-Stokes 
equation, but differs substantially at next level of the Reynolds' stress equation. 

8. Concluding remarks 

On the basis of the Boltzmann formalism, namely, nonequilibrium statistical mechanics 
designed for turbulence, two sets of equations are derived to comprise a closed set. The 
one is the group of equations governing Reynolds-averaged fluid quantities, and the other 
is the variable-separated version of the fluctuation-correlation equations that reduces to 
the classical Karman-Howarth equation in the homogeneous and isotropic limit. The 
two sets of equations are coupled through turbulent transports such as Reynolds' stress 
and turbulent heat-flux density. 

The key role of the first principle lies in that 

(i) for the Reynolds averaged equations, it gives a firmer basis than via phenomenologies, 
whereas, 

(ii) for the turbulent-transport equations, it reveals a hidden kinematical prerequisite 
for any turbulence governing equations to fulfill, thereby enabling to single out the 



17 



unique form. It is also shown that once their unique form has been established the 
whole sets of equations are able to be reconstructed a posteriori using phenomenologies 
alone. 

The latter set of equations is converted into a variable-separated form, leading 
to those governing solitary-wave function through which all the turbulent-transport 
properties are calculated. 

The equations are solved for a self-similar turbulent mixing layer, leading to a 
solitary-wave type solution in the physical-plus-eddy space. Reynolds stress is obtained 
through a simple integration of the solution over the eddy space in a form free from any 
empirical parameters, yet showing satisfactory agreement with existing experiment. 
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Figure captions 

Fig.l Solitary wave <7i(si,S3,0) on the plane of r\ = for two choices of the phase 
velocity c ; 

a) c = uo , where «o = [u(oo) + u(— oo)]/2 is the propagation velocity of eddies as 
observed by flow visualization, 

b) c = u (Taylor's hypothesis). 

Fig. 2 Reynolds' stress as calculated from Eq.(40) [18] using the solitary- wave solutions 
(Figs. la and lb), and compared with existing experiment [19]. 




c = u 



